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SUMMARY 


This report describes a theoretical investigation of three-dimensional 
relativistic klystron interaction. The relativistic axisymmetric equations 
of motion are derived from the time-dependent Lagrangian function for a 
charged particle in an electromagnetic field. An analytical expression of 
the fringing RF electric and magnetic fields within and in the vicinity of 
the interaction gap and the space-charge forces between axially and radially 
elastic deformable rings of charges are both included in the formulation. 
This makes an accurate computation of electron motion through the tunnel of 
the cavities and the drift tube spaces possible. The method of analysis is 
based on a Lagrangian formulation. Bunching is computed by using a disk 
model of an electron stream in which the electron stream is divided into 
axisymmetric disks of equal charge and each disk is assumed to consist of a 
number of concentric rings of equal charges. The individual representative 
groups of electrons are followed through the interaction gaps and drift tube 
spaces. 


INTRODUCTION 

The klystron is one of the most versatile electron devices used for 
amplification and generation of energy at microwave frequencies at high 
power levels. It has found many applications in communications, radar, mi- 
crowave energy sources for particle accelerations, microwave heating, and 
industrial processing. A recent advance in the depressed collector design 
(ref. 1) has made the high-power klystron a feasible high-efficiency micro- 
wave power source which is specially suited for transmission of large micro- 
wave power from space. 

Efforts to compute rigorously relativistic three-dimensional axisymme- 
tric electron motion in the klystron have been pursued at the Lewis Reseach 
Center. The present investigation serves as a continuation of such efforts 
with the objective of establishing a complete theory of relativistic kly- 
stron interaction. 

In this study, the Lagrangian formulation of a hydrodynamic beam model 
is used. The electron stream entering the interaction gap is subdivided 
into representative charge groups. The individual charge groups are fol- 
lowed through each interaction gap and drift tube space until the output 
interaction gap is reached. Thus electron overtaking and crossover are ap- 
propriately dealt with. The electron stream is divided into N axisymme- 
tric disks of equal charge per rf period, and each disk is assumed to con- 
sist of R concentric rings of equal charge. The rings are elastic and de- 
formable in the axial and radial directions, and the disks are assumed to be 
thin and the rings narrow. The velocity of each ring, its phase with re- 
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spect to the cavity gap voltage, and its radius at a specified position in 
various interaction gaps and in various drift tube spaces are computed from 
three-dimensional large signal programs. 

With the space-charge effects on the bunching process included, the 
debunching effects of the space-charge forces are appropriately dealt with; 
furthermore, with the radial motion considered, radial beam loading and the 
effects of a multidimensional charge motion on the field-electron interac- 
tion process are accurately formulated. 

This investigation separates conveniently into three parts: 

(1) Formulation of relativistic equations of motion in an electromag- 
netic field 

(2) A study of field-electron interaction in the first and second cav- 
ity gaps 

(3) A study of the large signal energy exchange process in the third 
and succeeding cavities up to and including the output cavity 


FORMULATION OF RELATIVISTIC EQUATIONS OF MOTION 


The relativistic equations of motion of a charged particle in 
linear coordinate system, such as the circular cylindrical system, 
safely derived from an invariant formulation of particle dynamics, 
the Hamiltonian or the Lagrangian. Mathematically these functions 
equivalent, but the Lagrangian method is somewhat more direct. In 
lows, the Lagrangian formulation is used. 

The Lagrangian of a particle of rest mass mQ and charge e 
electromagnetic field determined by the potentials i and A is 


a curvi- 
can be 
such as 
are 

what fol- 
in an 


L(t) = -hIqC^ y'l - 3^ - e-t + e(u • A) (1) 

where b = u/c, as usual in relativistic formulas. If q-j, i = 1,2,3, are 
the three coordinates of any generalized coordinate system, which define the 
position of the charged particle, then the Lagrangian dynamic equations of 
motion are 


d aL 




i = 1,2,3 


(2) 


where aL/aq-j = p-j are the canonical momentum components, which together 
form the total momentum vector of the particle. In circular cylindrical 
coordinates, 

_ .A 

u = ra^ + rcpa^ + za^ 


and 


u 


2 
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The field vectors in free space or in vacuum are given by the following re- 
lations: 


2 





B = V X A 


By performing the indicated differentiations, we obtain 

p = ^ = mr + eA 
3r 

p^ = mz + eA^ 

2 * 

p = mr <p + erA 

<p 9 


(3a) 

(3b) 

(3c) 


Of these equations, the first two, (3a) and (3b), are linear momenta, and 
the third, (3c), is an angular momentum. They are the sums of the relati- 
vistic mechanical momenta, differing from the ordinary (nonrelativistic) 
momenta only by the Fitzgerald factor which appears implicitly in the rela- 
tivistic mass m, and of a momentum of electromagnetic origin. It is as- 
sumed that the acceleration momentum is not important in the formulation. 
The relativistic mass m of the particle at velocity u is given by 



(see appendix A). Next, the derivatives aL/aq^ are given by 
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(4b) 
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In the case of axisymmetric field, the field is rotational ly symmetrical, L 
is independent of the azimuthal angle 9 , and we obtain from equation (4c) 

- 0 (4d) 

a 9 


The relativistic dynamic equation of motion of the charged particle can 
now be obtained by substituting equations (3) and (4) into equation (2) with 
the help of the following relations: 
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(5a) 
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B = 7 X A 
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When the result is regrouped. 


^ (rnr - mri^) = + rie^ - 

If mz = e^E^ + rB^ - r<pB^) 


(5b) 

(5c) 


We notice that expansion of the left side of both equations (5a) and (5b) 
will contain a term dm/dt; thus, when the energy conservation relation (eq. 
(B-4c) of appendix B) is used, equations (5a) and (5b) become (no = e/mo) 





+ 
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(6a) 
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(6b) 


Consider first the azimuthal equation of motion, equation (5c); the 
momentum integral yields 


P = constant 
9 


or 


If electrons have started from the cathode, where r = r^, 67= Bn, 
and 9 = = 0, the total amount of magnetic flux passing tnrougn a 

circle of radius rj in plane 1 (fig. 1) is found by integrating the 
product of the axial magnetic field B^ and the differential ring-shaped 
area 2irr dr from zero to r^ as follows: 
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r(v X A)^ dr 


where 62= (v x A)2 and in cylindrical coordinates is given by 

‘'•V 

for a/ 3 (p= 0 . We can now interpret the term r2Aj as the magnetic flux 
enclosed within the radius r^. Similarly the term is to be inter- 

preted as 172TT times the magnetic flux passing through plane 2 enclosed 
within a radius r2« Thus equation ( 7 ) can be rearranged to yield Busch's 
theorem 


or 


f \ 6^ (f y) (8) 

2irr'^ ^ 


where ’ll = '•'q is the cathode flux and equation (8) is now referred to any 
plane perpendicular to the axial direction. 




Figure 1. - Formulation of Busch's theorem. 
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Equations (6a), (6b), and (8) form a system of three-dimensional rela- 
tivistic equations of motion of an electron in an electromagnetic field, and 
they are rewritten as follows: 


r = 







(9a) 



c2 r 9 



(9b) 


9 = - 





(9c) 


where the E's and B's are the electric and magnetic field components 
which may exist in both the interaction gaps and the drift spaces in a mul- 
ticavity klystron. In an accurate treatment of klystron analysis, both the 
cavity circuit fields and the space-charge fields must be included in using 
equations (9). Thus 


E 

E 


r 


z 


= E + E 

r-cct r-sc 

(10a) 

^z-cct * ^z-sc 

(10b) 


where Er_gct E 2 _cct the RF cavity gap fields at the beam position 

in the radial and axial directions, respectively, and Ej._sg and E 2 _sc 
the radial and axial space-charge fields, respectively. 

Analytical expressions of the fringing RF electric and magnetic fields 
within and in the vicinity of the interaction gap obtained from a previous 
investigation (refs. 2 and 3) result in the following: 


and 


E^(r,z,t) = EqF^ exp(ju)t) 
E^(r,z,t) = EqF^ exp(jo)t) 






(lla) 

for -i. < z £ i 

(llb) 


E^(r,z,t) = EqG^ exp(jo)t)^ 
E^(r,z,t) = EqG^ exp(j(ot)J 


(12a) 

for 1 2 I 

(12b) 


where 
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r 




in which Eq is the electric field amplitude specified at the midplane of 
the cavity gap (z = 0 and r = a) and is related to the gap voltage by the 
following relations: 


V 


n 



cosh m(z 


fc)dz = 2 



sinh(m)i) 
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where is the gap voltage of the n^h cavity and m is the field-shape 
parameter. For a uniform field, m = 0. 

The three-dimensional nonrelativistic space-charge fields have been de- 
rived (ref. 3) in terms of two static Green space-charge potentials and 

S^p. The three-dimensional relativistic space-charge fields may be readily 
obtained by relating these Green function potentials as static potentials in a 
moving coordinate system and by applying a proper Lorentz transformation. 
Furthermore since_in the moving frame of reference the vector potential 
function A' or = 0, the fields in the laboratory frame can be found as 
follows: 
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where 
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are the two modified Green function potentials. In these potential expres- 
sions, 5 and p are the two normalized coordinates at the field point, and 
Co and pQ are the two normalized coordinates of the source points. 

The other symbols are defined as follows: 


Io(^n)» ^l(^n) 

^n 

k 

Pn 


modified Bessel functions of zero and first order, respec- 
tively 

Bessel functions of zero and first order, respectively 
nth root of Bessel function Jo(Xn) = 0 
free space wave number, u/c 

tunnel radius 


When the expressions for the cavity gap fields and the relativistic 
space-charge fields (eqs. (10) to (15)) are substituted into equations (9a) 
and (9b) and the resultant equations are normalized with respect to the 
tunnel radius a, we obtain the relativistic three-dimensional equations of 
motion in dimensionless form as follows: 

Normalized axial equation of motion: 



9 



Normalized radial equation of motion: 




r^ radius of cathode 

uq relativistic dc beam velocity, uqoRu 

Uqq nonrelativistic dc beam velocity, |/2(e/mQ)VQ = 5.93x10^ m/sec 

9 5 

Veq equivalent beam voltage, m^c /e = 5.11x10 volts 


Be “/uqo 

<p 0)t 

I normalized axial coordinate, z/a 

i dC/de 

p normalized radial coordinate, p 

p dp/de 
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'?(p.c) 



p dp 


Vq cathode flux, ira^Bc 

tocO nonrelativistic cyclotron angular frequency, (e/mo)Bo 


y Wmniz 

n — ~ 

00 

We note that equations (16) and (17) are highly nonlinear, cannot be solved 
analytically, and must, therefore, be solved on a high-speed digital compu- 
ter. Furthermore the modulation index a should be specified before equa- 
tions (16) and (17) can be used. 

For the input cavity gap (the first) and the second cavity gap a is, 
by definition, small (e.g., <0.02 and <0.07, respectively). In particular, 
for the first cavity, ai can be related to the power input to the cavity 
and to the total shunt resistance of the cavity by the 

relation 


a 



f- 


where Rgh = Rcav Rb RexJ Rcav is the shunt resistance component due to 
cavity losses; R^ is the component from the beam loading, which may be 
negative but is ordinarily positive and often the dominant component of the 
total shunt resistance; and Rg^ is the reflex load resistance due to 
external losses (load) coupled to the cavity (zero in the case of the input 
cavity). Generally Rgh is related to the unloaded Qy ana the R/Q 
value of the cavity by the relation 


n 

*^u = R7Q 

Hence the value of Rgh can be obtained by values of and R/Q 
which are usually the design parameters for a given cavity. 

CALCULATION OF BUNCHING CURRENT AND VELOCITY MODULATION IN SECOND CAVITY 

Interaction of an unmodulated electron stream with the input cavity 
(first cavity) gap fields gives rise to the velocity modulation of the 
stream. As the beam moves along the first drift tube, the velocity varia- 
tion of the stream is converted into density modulation, and the electron 
current at any point in the drift tube space can be computed from the kine- 
matic of the electrons. 


11 



ElECTRON BEAM 



If we divide the beam into a number N of disks of electrons and we 
divide each disk further into R elementary charge rings (fig. 2), an 
elementary charge ring 2irro dro dzQ containing a charge density 
oo(zo,ro,to) at time tQ at a later time t becomes p(z,r,t) 
in an element of charge ring 2irr dr dz. From the charge conservation 
principle, we relate these charges by the following equation: 

p(z,r,t)2irr dr dz = Po(2o’'^0’^0^ EF”'o '"'o 

NR 

The modulus is taken, since the sign of the charge is unchanged even if the 
electrons overtake one another. The summation over (2irrQ dro dzQ) implies 
that all elements of the charge ring which entered the first cavity gap ap- 
pear within the charge ring 2irr dr dz at the later time t as a result of 
overtaking and trajectory crossing. Thus 


p(z,r,t) 





•^0 '*'^0 ^^0 
r dr dz 


(18) 


where the charge density po(zo,ro,to) is assumed constant over the beam 
cross section at the injection plane z = Zq, and, in terms of the beam cur- 
rent Iq, 


I 


0 


■ f,2„ 

Tfb u. 


We note that rg is the mean radius of the charge ring at zq. The linear 
charge density in any annular ring and at any given displacement plane 
transverse to the direction of the beam is written as 
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da(z,r,t) = -2i;p(z,r,t)r dr 


_ -2trPQ(ZQ,l"Q»tQ) 


^ |’'o ‘"'o^ 


NR 


by equation (18). 

Furthermore, by expressing 



as 


^^0 



zr 

“u^XzTrvt) 

dt 


u^ir,r,tj “ 


-wr 


(19) 


we can write equation (19) as 


da(z,r,t) 


-2irpo(Zo,ro»to) 


u^U.r.t) 




(20) 


In the computation of bunching current in the second cavity, it is assumed 
that there is no net current in the radial direction (i.e., small signal 
theory is assumed); then the incremental current through the annular ring 
can be written as 


i(z,t,t) = -27Tp(zQ,rQ,tQ)u2(zQ,rQ,tQ) f 

Jo 
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( 21 ) 


In terms of the normalized parameters l and p, equation (21) can be 
written as 


i( 5,P, 0) = 21 




b/a 


ZP 




( 22 ) 
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It is to be noted here that pq = ro/a is not the dc beam charge density 
used prior to equation (22). 

The total electron convection current may be conveniently expressed in 
terms of a Fourier series written in the variables 0 as 


in(C,P,e) = Aq + ^ (A^ cos ne + sin ne) 

n=l 


(23) 


where 



(24) 
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n 



■'p( e)cos n 0 d e 




cos no d 0 Q 


(25) 


and 


B 


n 



■'p( 5,p, 0)sin no do 


sin ne de^ 


( 26 ) 


in which o is the phase of the electron at any point along the stream. 
This is to be related to the injection phase of the electron at the 
time of entering into the interaction gap space of the first cavity. The 
interaction gap space is defined (fig. 3) as extending from z = -( a + 2a) 


14 



to 2 = +( n + 2 a), where 2 a is the tunnel diameter and 2i is the 
cavity physical gap length. Thus 


t - t * iii-LiSi * 

\ - ^0 u„ 


( i + 2 a) 


or 





®9 


uL^ - ( A + 2 a) 


(27a) 


In terms of the normalized parameter K, this can be written as 

®2- »g*(a-r-2) 


(27b) 


where U 2 I is the electron velocity at the plane of exit of the first 
cavity gap to be either computed or substituted from small signal theory. 

The bunching current at the midplane of the second cavity gap (z = 22 ) 
can now be obtained by, first, letting e = 03 and, second, substitut- 
ing equation (27b) for e in equations (24) and (25). The result is 



CAVITY GAP 1 CAVITY GAP2 CAVITY GAP n 



Figure 3. - Geometrical parameters of multicavity klystron. 
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The ratio of harmonic current amplitudes to direct current is given by 



for n = 1, 2, 3, . . . . The phase angles of the n^h harmonic with 
reference to a pure cosine or pure sine wave of the same frequency are given 
by 


A 


n 



(30) 


Once the electron convection current is obtained, the induced current, 
and hence the induced voltage, can be computed. Since 02 is small in the 
second cavity, the current induced there as a consequence of the bunched 
electron beam traversing the cavity gap can be computed by the following re- 
lation: 


^(ind)2 


M^M^i2 ( C»P j ^2 ) 


(31) 


where i2(5jPf02) the electron convection current passing through the 
midplane of the second cavity, and and are radial and axial 
beam coupling coefficients, respectively, given by the following relations 
(ref. 4): 


and 


M = 
r 




TqIV) 


sin 3 2 i 

-371- 


(31a) 


(31b) 


where 



is the radial propagation constant, = u/c. 
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Be = (o/UQ, and Iq and are the modified Bessel functions of the zero 
and first order, respectively. We note that the axial beam coupling coeffi- 
cient is to be evaluated at the edge of the cavity gap. The voltage across 

the second cavity gap (i.e., across the capacitor of the equivalent shunt 
circuit, fig. 4) is given by 


'^gZ " ^2^(ind)2 

where Z 2 is the dynamic shunt impedance of the cavity given by 

^2^“^ " 1 + 2j6Q^j 


in which Rgh is the shunt resistance at resonance, and s is the frac- 
tional deviation from resonance and is defined by 
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(32) 


where the subscripts 2 imply that these quantities are to be evaluated in 
cavity gap 2. 
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Figure 4. * Schematic of multicavity klystron. 
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LARGE SIGNAL ENERGY EXCHANGE AND CALCULATION OF INDUCED CURRENT 

In a multicavity high-power klystron, interaction between the bunched 
electron stream and the cavity gap fields in the third and succeeding cavi- 
ties up to the output cavity becomes complicated. The simpler method of 
computing the induced current used in the section CALCULATION OF BUNCHING 
CURRENT AND VELOCITY MODULATION IN SECOND CAVITY is no longer valid. In 
this section, a general method of approach is used. This method is based on 
Shockley (ref. 5) and Ramo (ref. 6) theory. 

Consider an electron of charge -e that travels a distance dl inside 
the cavity interaction gap space (fig. 3). The gain or loss of its kinetic 
energy is given by 


dW = -e 



dl 


(33) 


By the theory of Shockley and Ramo, the exchange of energy is related to the 
induced current that flows in the external circuit (i.e., the cavity) and 
the induced voltage that appears across the cavity gap by the following re- 
lation: 


dW = -I 


V 

ind ind 


dt 


(34) 


Setting equation (33) equal to (34) yields 


I 


ind 





(35) 


where Ecct is the cavity gap field (the circuit field) and u = dl/dt is 
the electron velocity. 

Since many electrons are present inside the cavity interaction gap 
space at any instant of time and the electron transit time within the inter- 
action gap space is finite, each elementary current carried at a given mom- 
ent of time by each elementary charge within the interaction gap space will 
contribute to the total induced current. These elementary currents in the 
plane over which the beam enters the interaction gap space, and in the plane 
located some distance away in the direction in which the beam is traveling, 
differ not only in velocity but also in phase. By summing up the elementary 
current or charge within the interaction gap space, we obtain the total in- 
duced current 


J ind 

T 

where p(r,z,e,T) is the charge density inside the interaction gap space. 
Alternatively, in terms of convection current density Jc (= p^» equa- 
tion (37) can be written as 
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(36) 
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ind 


(t) 




dx 


where dt = r dr d<p dz is the volume element. The integration is taken 
over the whole volume of space of the interaction gap space occupied by the 
bunched electrons at the time t. We note that, since the cavity gap field 

Ecct proportional to the induced voltage V-jpH. for a given electron 
velocity u, the induced current is independent of the voltage across the 
cavity gap. 

Using the first expression, equation (36), we find 
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i nd 


(t) 



b'(2) 


p(r,z,t)u 


cct 

ind 


r dr dz 


(37) 


where b'(z) is the edge of the beam as a function of z 


and is defined by 


b'(z) 



(38) 


in which rg is the radial position of the electron at the entrance to the 
interaction gap space, here to be referred to as the entrance plane of the 
first cavity gap space. 

The integral of equation (37) can be evaluated by applying the charge 
conservation principle (see fig. 2) to write the charge density p(r,z,t) as 
function of the phase of the injected electrons at the entrance to the in- 
teraction gap space of the first cavity as follows: 


p(r,z,t) 





'""o 


“"o 

r dr 


dz ' 






r dr 


u^(r,z,t)dt 




*^eq 


- d6Q 





(39) 


where beq is the equilibrium, or the average, radius of the outer charge 
ring centroid given by the relation 
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^/2 




T72 


,1/2 


(1 - 3^) 


“cO 


in which wpo and ucO ^re the nonreiativistic plasma and cyclotron 
angular frequencies, respectively, and R^, is the relativistic reduction 
factor as given in the section FORMULATION OF RELATIVISTIC EQUATION OF 
MOTION. 

Substituting equation (39) into equation (37), we obtain 


(a+2a) -b'(z) 




M.I 


- ( 1+2 a ) 0 


u^(r, 2 ,t) 


^cct 

Vind 


= 2irl, 


eq 
b' (z) 



(ii+2a) -b'(z) 


/ 


-{i,+2a) 0 


^r-cct 

'^ind 



(40) 


where U 2 (ro,zo,lo) = '^0 velocity. In terms of the normal- 

ized parameters c and p, this is written as 




b' (z) 




. 2 . 1 , 


Mf f 

^ U+M n 


"r-cct 

''icid 


+ 


^z-cct ] 

^•nd / 


X 



a^p 


dp dz 


(41) 


Furthermore the induced current I-jpcj may conveniently be represented by a 
Fourier series in the variable e as follows: 
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n=l 


cos ne 

IT 



ne de 
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lind 


(42) 


where lind(t) is given by equation (41). 

The two integrals in equation (42) can be evaluated, first, by trans- 
forming the variable e in terms of eo» electron injection phase at 
the entrance plane to the interaction gap space of the first cavity. With 
reference to figure 3, we note that the times, tj, t?, t 3 , . . ., t^, taken 
by electrons to pass through the rf gaps are interrelated: 


t - t. + - 

3 2 u_i u 


zl 


z 2 


-t . 2 (»* 2a) ^ '-1 ^ '•2 ~ 

0 “O “zl "z2 


2a) 


In terms of the normalized parameter 5 , this can be written as 



In general 


n =2 



where 0g = 2(a+ 2a)(w/uo) is the dc transit angle of the electron in tra- 
versing the first interaction gap space. 

Next, by relating the cavity gap field Tcct with the induced voltage 
Vind through equations (11) to (13), we can write 
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""'ll ”"<"n ■" "n) 


2 s1nh(ml) Vp cos(0^ + 


m 


2 sinh(ms,) yF 

Furthermore, for simplicity, equation (42) is written in terms of the Four- 
ier coefficients, Aq, Ap, and Bp. By a substitution of equation (43) 
into equation (42), we obtain 


An = 



[?f^l 

2 s1nh(m*,) 




(fc+a) ^ b'(z) 


X p dp d5 d0Q (44) 


An = 


ml^a^ 

r ~l 

eq 

2 sinhOiiji) 

b'(z) 


(fc+a) . b'(z) 


Bn = 


mloa^ 

^eq 

2 sinh(ma) 



n -"T" r 

d. 

Jo J.i^Jo 


X cos p dp dc dSn (45) 


X sin p dp dc de_ (46) 


where 



X = -e- 

The G's are given bv equations (11) to (12), and the variable e is 
given by equation (43). The amplitudes of the harmonic to dc can be ob- 
tained by using the relation 


(47) 


^(ind)n 
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A + B 
n n 


(48) 


and the phase angle is related by 


+ -1 


(49) 
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With the induced current obtained from equations (42) to (49), the mod- 
ulation index for the third and succeeding cavities up to the penultimate 
cavity can be computed through an iteration procedure. To illustrate this 
method of approach, let us consider a case of computing 03 (i.e., the gap 
voltage of the third cavity) as an example: 

(1) Find the induced current as computed at the midplane of the third 
cavity gap, called I(ind) 3 - 

(2) Let the gap voltage V „3 = 0, and compute the gap voltage using the 
relation 

'^g3a ~ ^3^“^^(ind)3 


where 


Z3U) 


R 


sn 


1 + 2j6Q 


u3 


(3) Using the value of Vg 3 g as obtained in step ( 2 ), formulate and 
solve the equations of motion, equations (16) and (17), and thus obtain the 
induced current, called 139 . 

(4) Using i 3 g, as in step (2), compute the gap voltage, called the 
gap voltage Vg 3 b. 

(5) Repeat steps (3) and (4) until a gap voltage Vg 3 |^ is obtained so 
that the following condition converges: 



In the case of the output cavity, the gap voltage VgQ is simply 

V W"»( 1 nd )0 <5<» 


where l(ind )0 "is the current induced in the output cavity, and Zout(‘^) ‘is 
the dynamic shunt impedance of the output cavity given by 


■7 ! \ ^shO 


(51) 


Power delivered to the external load of the cavity is obtained by 

1 


out 


2 '^®^ind)o''gO 


(52) 


where Vgo is the complex conjugate of VgQ. Finally, the output 
cavity efficiency can be computed in the usual way: 
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„ Pout (l/^)P^-nnd)o'' 5 o ,53 

" ^ ''o'o ‘ 

where Vq is the dc beam voltage, and Iq is the dc beam current, which is 
related to the perveance by the relation 


K 


r 




( 54 ) 


where Kn = Iq/Vq^^ is the nonrelativistic perveance, and Vg = moc^/e 
= S.llOxIO^ volts is the equivalent beam voltage. 

SUMMARY 

Three-dimensional equations of motion of relativistic electrons in an 
axisymmetric field in the presence of space-charge forces has been ob- 
tained. The field-electron interactions in various interaction gaps and in 
various drift spaces are investigated by using a Lagrangian formulation of 
the hydrodynamic model. The equations of motion are used to find the velo- 
city of each disk of electrons, its phase with respect to the interaction 
gap voltage, and its radius at a specified position in various interaction 
gaps and in various drift spaces. The velocity modulation and electron 
bunching in the input and second cavities are formulated by using the small 
signal approach; however, both the radial and axial coupling coefficients 
are useo to compute current. The induced current in the third and succeed- 
ing cavities up to and including the output cavity is computed on the basis 
of Shockley and Ramo theory. Power and efficiency are calculated in the 
usual way. 
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APPENDIX A 


MODIFICATION OF PARAMETERS DUE TO RELATIVISTIC VELOCITY 


Mass Increase 


When the electrons are accelerated through very high potentials, they 
experience an increase in mass which modifies the velocity-potential rela- 
tionship. The increase of the electronic mass with increasing electron 
velocity is given from the transformation law originally proposed by H. A. 
Lorentz (1904): 


m = 



(Al) 


where e = u/c, as usual in relativistic formulas; c = 3x10® meters per sec- 
ond is the velocity of light; m is the mass of the electron at the velocity 
u; and mo is the electronic rest mass at zero velocity. 

Velocity of Electron 

If W is the total energy of the electron in joules, the relation 


W = 


me 


(A2) 


expresses the equivalence of the instantaneous mass and the inertia of 
energy. The total energy of an electron is the sum of the energy inherent 
in its mass plus its kinetic energy. The total energy change of an electron 
when it is accelerated through a potential of Vq volts, can, according 
to equation (A-2), be expressed by 

2 

eVg = aW = c (m - m^) 




0 


(A3) 


Thus 



V„ + 1 


mQC 


2 '0 


(A4) 


By introducing the energy equivalent of the electron rest mass mo as 


eV 


e 



(A5) 


we obtain 
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e 



0.5110x10^ volts 


(A6) 


Hence equation (A4) can be written as 


m 




(A7) 


When this is substituted into equation (Al), we obtain the relation between 
the electron velocity u and the accelerating potential V: 


u 


0 " 



(A8a) 


or, alternatively. 



. W(V + 1022) 

^ V + 511 


= UoqRu V in kilovolts) 


(A8b) 


where 

“00- V^^''o 


is the nonrelativistic expression of the electron velocity, and 



is the relativistic correction factor. 

Angular and Plasma Frequency of Electron Stream 

The angular frequency of revolution of the electron in the homogeneous 
magnetic field B is defined as the cyclotron frequency given by 
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At moderate electron velocities, the increase in electronic mass is negli- 
gibly small; hence a»c is considered to be independent of the velocity of 
the electron. However, this is not true for the high-velocity electron. 

The cyclotron frequency for large velocities follows by substitution of the 
relativistic mass of equation (Al) for m; 


(0 


c 






(AlO) 


where wco = (e/mo)B is the nonrelativistic cyclotron frequency. In 
terms of Vg, equation (AlO) can be written as 


0 ) 




(All) 


In like manner, the relativistic plasma frequency can be found by substitu- 
ting the relativistic mass, equation (Al), and the relativistic velocity, 
equation (A8b), for m and u, respectively: 



CoTTUq 


= (0 


po 



- 1/4 


where a> « = 1 / r , — "is the nonrelativistic plasma frequency, and u^n* as 

pu y SqttUqq uu 

given in equation (A8b), is the nonrelativistic velocity of the electron. 


27 



APPENDIX B 


ENERGY CONSERVATION PRINCIPLE 

Consider an electron which has been accelerated through a potential 
Vq volts and gained kinetic energy in the amount of oVq; then, by the 
equivalence of mass and energy principle, the change in the total energy of 
the charged particle is 


fiW = eV 


0 


2 2 
c (m - m„) = c 6m 


(Bl) 


by equation (A3) of appendix A. 

Work ^ne by the force in moving the charged particle of mass m a 
distance dr is given by 


dW = 7 • dr = e(T+ lT X ?) • ^ (B2) 

where the Lorentz force law has been involved on the right side of the 
equation. Now, power is the rate of change of energy; thus, dividing both 
sides of equation (B2) by dt, we obtain 

^.e(rw,B) 4 


= e(? • T) (B3) 

where IT = dr/dt, and ^7x7)* IT is identically equal to zero. Equation 
(Bl) is differentiated with respect to time, and the resultant expression is 
used to write equation (B3) as follows: 


or 



(B4a) 


(B4b) 
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In a cylindrical coordinate system, this is written as 






% 


d 

dt 




<P 



(B4c) 
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